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Abstract. The results of Coulomb gauge and Landau gauge lattice QCD 
simulation do not agree completely with continuum theory. There are indica- 
tions that the ghost propagator in the infrared region has strong fluctuation 
whose modulus is compatible with that of the color diagonal ghost prop- 
agator. After presenting lattice simulation of configurations produced with 
Kogut-Susskind fermion (MILC collaboration) and those with domain wall 
fermion (RBC/UKQCD collaboration), I investigate in triple gluon vertex 
and the ghost-gluon-ghost vertex how the square of the color antisymmet- 
ric ghost contributes. Then the effect of the vertex correction to the gluon 
propagator and the ghost propagator is investigated. 

Recent Dyson-Schwinger equation analysis suggests the ghost dressing 
function G{0) = finite and no infrared enhancement or ac = 0. But the ghost 
propagator renormalized by the loop containing a product of color antisym- 
metric ghost is expected to behave as < cc >r= — '^^^ ^ with G{q^) oc g~^"<3 
with ac = 0.5, if the fixed point scenario is valid. I interpret the ac = solu- 
tion should contain a vertex correction. The infrared exponent of our lattice 
Landau gauge gluon propagator of the RBC/UKQCD is ao = —0.5 and that 
of MILC is about -0.7. 

A possible interpretation of the origin of the fiuctuation is given. 



1 Introduction 

Color confinement and chiral symmetry breaking are the most important char- 
acteristics of the infrared (IR) QCD. We studied their mechanism by measuring 
the gluon propagator DA{q'^), the ghost propagator Dciq^) [H El El HI [U [6] 
and the quark propagator [7J in lattice QCD and comparing results with those 
of Dyson-Schwinger equation (DSE). From the condition on the IR fixed point 
of the running coupling in the Landau gauge which is calculated as as{q'^) = 
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q^Dciq^^^DAiq^)-, the ghost dressing function defined as G{q^) = q^Dciq^) 
in d— dimensional system was shown to behave as G{q^) oc [8j and 

q^DA{q^) oc g'^^^ ^ith 

A-d 

2aG + aD + ^— =0 (1) 

in the IR, or it is more singular than the free case. 

Recently an SU(2) gluon propagator of large lattice was shown to be in- 
frared finite[9] which was in conflict with the prediction of the DSE[101lllj which 
showed that it is infrared vanishing, or the IR exponents which is defined as 
q^DA{q^) oc q~'^°'o \^ smaller than -1. The results could be attributed to the finite 
size effect, since infrared vanishing gluon propagator was predicted also from the 
continuum theory of Zwanziger [121 [T3l [HI |T3] which considered uniqueness of 
the gauge field or the celebrated Gribov problem |15j and proposed restriction of 
gauge configuration to the fundamental modular region, and predicted without 
incorporating two-loop contribution k = uq = 0.59. 

In the DSE approach, by incorporating the two-loop gluon contribution, an 
IR finite (0/5 = —1) gluon propagator was also proposed [16j, although the 
incorporation of the two-loop or squint diagram in this work is not without 
ambiguity [T7]. Recently a new solution of DSE with IR exponent = — 1 for 
the gluon propagator but ac = for the ghost propagator was proposed p^[T9]. 
In this case, the running coupling in the Landau gauge vanishes at g = 0. This 
behavior was observed in our lattice simulations [3] but we concluded that it is 
an artefact. There are arguments against the IR exponents of the new DSE |20| . 

One should ask the validity of the definition of the running coupling which 
is based on the tree approximation. I observed in \21] that the loop corrections 
through the color antisymmetric ghost could affect IR features of QCD. The 
color antisymmetric ghost and its relation to ghost condensates was discussed in 
[22] and the upper limit of the modulus of color antisymmetric ghost is measured 
in SU(2) [231 El] and in SU(3) [3 [6]. 

The running coupling extracted by the JLab group from the experiments 
shows freezing to a value close to 3.14 |25j . and the conformal field theory based 
on the Crewther relation predicts 0^(0) = vr [26^ I27|. A lattice simulation of un- 
quenched configurations in Coulomb gauge is consistent with the JLab extraction 
[28]. 

In this paper I summarize our lattice simulation data of the ghost propagator 
in sect. 2 and discuss the effect of color antisymmetric ghost in the gluon- ghost- 
ghost vertex in sect 3. Since the color antisymmetric ghost could modify the 
ghost-gluon vertex, I discuss a reconsideration of the Slavnov identity due to the 
presence of the color antisymmetric ghost in sect 4. A discussion and conclusion 
are given in sect 5. 

2 Ghost propagators in lattice simulations 

In Landau gauge {d^A^ = 0) and in Coulomb gauge {diAi = 0), we adopt the 
logC/ type gauge field A, i.e. U^,^ = e^^^'^*, aI^^ = -A^^^, and for the gauge 
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Fu{g) = W= tr(^V^|,J. (2) 

x,/i=l, 2,3,4 

in Landau gauge and 

x,i=l,2,3 

in the Coulomb gauge. 

The gauge field ^0(2^) in Coulomb gauge can be fixed by the following mini- 
mizing function of ^(xo). 



Fu[g] = \\A%f = Y^T(A^jAl,), (4) 



but in this work I leave the remnant gauge unfixed. 
In Landau gauge, the ghost propagator is defined as 

FT[D^^{x,y)] = FT{tr{A'^{{M[U])-'UA'), 

where A4 = —d^Dfj^. We solve the equation with plane wave sources. 

-d^D^f',{x) = ^A'smqx (5) 



d^D^f,{x) = ^A''cosqx. (6) 



The color diagonal ghost propagator is defined as 
Dciq) ^ ^ 



N^-IV 

xd''\{A''cosqx\f',{x)) + {A''smqx\f',{x))), (7) 
and the color antisymmetric ghost propagator is defined as 

xr''%{A''cosqx\f,ix)) - (A'^sinqxlf^ix))) 

(8) 

where = 2 for SU(2) and 6 for SU(3). 

In our lattice simulation we use the Kogut-Susskind (KS) fermion of MILC 
collaboration [29l [3^ and the domain wall fermion (DWF) of RBC/UKQCD 
collaboration [31], whose lattice size and parameters are shown in TABLE [TJ 

In [SI [6], we showed that the modulus of the color antisymmetric ghost prop- 
agator of quenched configuration is small and its variation is large and in [7| the 
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Table 1. The parameters of the lattice configurations. 
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Figure 1. The color antisymmetric ghost dressing function logj^Q(|0(g)|g'^) (Orange stars) and 
color diagonal ghost dressing function logG(g) (Blue diamonds) of quenched SU(2). (67 sam- 
ples) Slope defined by the lowest six points of the former is -1.8(1) and that of the latter is 
-0.45(2). 



color diagonal ghost propagator is essentially temperature independent while the 
color antisymmetric ghost is temperature dependent. 

In FigsH] and [2] we show the color diagonal ghost dressing function G{q^) 
and the color antisymmetric ghost propagator (j){q) multiplied by of quenched 
SU(2) after parallel tempering gauge fixing [T] and those of unquenched MILC3J 
configurations, respectively. Errorbars are standard deviations. 

The slope of the color diagonal and that of the modulus of the color anti- 
symmetric ghost depend on the number of flavors. The extrapolated fitting lines 
of SU(2) and MILC3/ cross at around q = 60MeV and MILC/(/3imp = 7.09) 
at around 30MeV (Figl2]). In the case of DWFqi, they cross at around lOOMeV 
(FiglSj) . Since the lattice size L of DWFqi is less than 2 fm, finite size effect on 
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Figure 2. The color antisymmetric ghost dressing function logj^Q(|0(g) |g^) (Red stars) and 
color diagonal ghost dressing function logG(g) (Green diamonds) of MILC/(/3imp ~ 7.09). (21 
samples) Slope defined by the lowest three points of the former is -1.4(3) and that of the latter 
is -0.27(8). 
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Figure 3. The color antisymmetric ghost dressing function logj^Q(|(/>(g) |(?^) (Orange stars) and 
color diagonal ghost dressing function logG(g) (Magenta diamonds) of DWFoi. (50 samples) 
Slope defined by the lowest three points of the former is -1.4(1) and that of the latter is -0.38(8). 



the color antisymmetric ghost |24j could be important in this case. If there are 
ghost condensates, the color antisymmetric ghost propagator is expected to tend 
to a constant in the IR. 

In the UV, the color diagonal ghost dressing function remains to be a con- 
stant, i.e. ac = even when the color antisymmetric ghost exists. In |2l|, I 
showed that color antisymmetric ghosts yield a dominant component of the prop- 
agator const/g^ in UV when the IR exponent has ac = 0.5. In the next section 
I discuss the gluon-ghost-ghost vertex. 

3 Color antisymmetric ghost in the gluon-ghost-ghost vertex 

In [21], I showed that in the triple gluon vertex and in the gluon-ghost-ghost 
vertex, product of color antisymmetric ghost produces real matrix elements when 
the product makes a color index in Cartan subalgebra. The vertex can be inserted 
in the ghost propagator as in FigjU or in the gluon propagator as in FigjSl 
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Figure 4. The dressing of the ghost propagator by the gluon. The dashed line represents a 
ghost, the thin line a gluon, and the thick line a quark. 




Figure 5. The ghost loop contribution in the gluon propagator. The dashed line represents a 
ghost and the thick line is a quark. 



In the ghost propagator, color indices h and d in FigU] can make a color 
antisymmetric ghost when the index a is in Cartan subalgebra, and thus the 
ghost propagator is not necessarily color diagonal and fluctuates by the operator 

On the contrary, in the gluon propagator, color indices e, d, h shown in Figl5] 
can be rotated. When the color index a is specified in Cartan subalgebra, color 
antisymmetric ghost can be chosen not only at d and h but also at d and e or h 
and e. Relative sign of the latter two cases and the former is random [21]. Thus 
the gluon propagator is effectively color diagonal. 

In lattice simulation, the Kugo-Ojima criterion is satisfied in the unquenched 
simulation but by about 80% in the quenched simulation. We investigated dif- 
ferences of the color SU(3) ghost propagator in Landau gauge and observed that 
the color antisymmetric ghost propagator in quenched configuration of 56'* lat- 
tice is random as shown in FiglH but that in the unquenched configuration its 
randomness is lost and its slope of the modulus as the function of momentum is 
steeper than that of the color diagonal ghost propagator and its extrapolation 
to IR becomes compatible with that of the color diagonal component. 

The product of the color antisymmetric ghost affects the triple gluon or 
ghost-gluon-ghost vertices and gluon propagator and the ghost propagator via 
loop corrections. The QCD effective coupling is affected by the ghost propagator 
and there appear differences in Coulomb gauge and in Landau gauge. 

In [21], I found that the ghost propagator with a gluon-ghost-ghost loop 
with a tadpole vertex between the two ghost propagators introduces an effective 
propagator of a const/p^ if the IR exponent of the ghost propagator ac = 0.5. 
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Figure 6. Scatter plot of the color antisymmetric ghost dressing function logj^Q(|(;/)(g)|g^) of 
quenched 5(7(3). (10 samples) 



It means that in the UV, the dominant component of the ghost propagator is 
const /p^ as observed recently by the Dyson-Schwinger approach [18l[19]. In the 
Dyson- Schwinger approach of [11] the exponent in the UV is the same as that 
in the IR. Our calculation of the loop diagram and the lattice simulation of the 
ghost propagator suggest that the exponent in the UV is but that in the IR 
is finite and enhances the singularity. If ac = 0, additional vertex singularity of 



-4k 



q yields the same const/(7 dominant behavior in UV. 



4 A comparison with continuum theory of infrared QCD 

The origin of the fluctuation of the ghost propagator could be the Gribov copy 
effect and I would like to discuss applicability of the Faddeev-Popov quantization 
in the infrared region. 

In the Faddeev-Popov quantization method, Zavialov showed in his 
book [32] that for an arbitrary functional F{A, c, c) depending on fields 
j4^(xi), c"(yj), c^{zf:) and let 5F{A, c, c) denotes its variation under BRST trans- 
formation, the Green function of 6F is zero, i.e. {S5F) = 0. 

Zavialov took F = A1^{x)c'^{y) and calculated {S6F) = {6FS). He showed 
that 



F{-inf'''A,{x)) 



SF 



6A-{x)' 



(9) 



which implies that Afj_{x) can be replaced by the functional derivative, and would 
be valid since gluon propagator is color diagonal. 

In the case of ghost, however it is not evident that the ghost field can be 
replaced by the functional derivative as 



F(-inc"(x)) 



SF 

6c'' (x) 



(10) 



Since gauge fixing can be defined only locally, gauge potentials are partitioned 
into patches corresponding to copies. The Faddeev-Popov quantization in the 
high energy region could be inefficient in infrared due to overlapping of patches. 
There are proposals of stochastic quantization [Sj [331 (SI [35], but whether the 
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Figure 7. The cancelling mechanism of longitudinal gluon and ghosts. 



unique gauge is attained through the stochastic quantization is not evident due 
to large fluctuation from color antisymmetric ghost. 

The argument does not mean the Slavnov identity is broken, since in the 
tree level, the color antisymmetric contribution cancels out among themselves. 
However, in the one loop, products of the color antisymmetric ghosts affect IR 
features of QCD. 

Similar argument applies to the renormalization of the gauge theory. In 1971 
t'Hooft showed that in the massless Yang-Mills field theory, the gauge invariance 
can be restored by incorporating finite number of counter terms including ghosts 
and longitudinal gluons in the system |36j . He showed, if the longitudinal gluon 
forms a pair of color diagonal ghosts and the loop intersects with the circle 
that denotes a set of particles on mass-shell, the contribution is cancelled by 
the cutoff denoted by A (Fig. 7). But when a pair of color antisymmetric ghosts 
are produced and there are Gribov copies, whether one can define A globally 
is not clear. In compact lattices, we observe fluctuation of propagators is large 
for momenta which deviate from the diagonal of the four dimensional system 
(cylinder cut region), which may be related to this issue. 

In a semiperturbative DSE calculation [37], and in a lattice simulation [38j , 
the infrared ghost-gluon vertex in Landau gauge was claimed to be close to that 
of tree level. These analyses are done in quenched approximation, where modulus 
of color antisymmetric ghost is random as shown in FiglH It is not clear whether, 
in unquenched systems, a product of color antisymmetric ghosts that appear in 
one loop level does not affect the theory. 

4.1 Kugo-Ojima color confinement criterion 

In 1971, Taylor [39j pointed out that the Ward identity in QED i.e. Z1/Z2 = 1 
can be extended to QCD as Z1/Z3 = Z^^j^/Z^. In 1979 Kugo and Ojima [iO] 
showed that if color symmetry is not broken Zi/Z^, = Zi/Z^ = Z^^j^/Z^. The 
IR exponent of the ghost changes the color confinement criterion of Kugo and 
Ojima, which says that in the Landau gauge a coefficient u°'^{(f') = 5°'^u{p'^) and 
if at momentum zero u(0) = — 1 the color confinement occurs. 
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When we define 

(cc) 



{{A,, X c)c)ip/ = -iq^J'iq'^) 

and 

{D^cc) = {d^cc) + {{A^ X c)c) = iq^l + ^(g2))_^_^ 
we obtain, provided c in Afj_ x c and c in x c do not couple color diagonal, 

{D^c{Af, X c)) = (5^cc)(c(^^ X c))ip/ + ((^^ X c){Au x c))ip/ 

/ (V.-^)(c(^mXc))ip/. (12) 

Although our finite size lattice simulation suggests ^(0) = 0, there could be 
a contribution of the ghost propagator of the type appearing between c and c 
that is not proportional to 5"''. Thus, l + n(0) = 1 + ^(0) is not necessarily equal 
to and ^(0) can deviate from -1. 

In the continuum limit, a solution with ^(0) =finite and ac = is possible 
Then 1 + ^(0) is not necessarily of higher order 0, as required by Kugo 



4.2 QCD effective coupling 

In Landau gauge, the running coupling as{q'^) can be calculated from the ghost- 
gluon coupling, the triple gluon coupling or the quark-gluon couping. The ghost- 
gluon coupling is given by the product of the ghost dressing function squared 
times the gluon dressing function. 

q^DciqfDAiq) oc as{q^). (13) 

The small ghost dressing function causes IR suppression of the running coupling, 
as shown in FiglHl 

In Coulomb gauge, the corresponding coupling constant aj{q'^) is obtained by 
choosing the interpolation gauge parametrized by rj between the Landau gauge 
{rj = 1) and the Coulomb gauge (ry = 0). Using the limit of an integration over the 

dqo ri\q\^ 3 



4th component of the momentum lim / 7i — ttttt = — , we obtain 

q'DG{qfD'XiQ)^M<l'), (14) 

where D^J{q) is the 3-dimensional gluon propagator |42j. No IR suppression oc- 
curs in aj{q), as shown in FigEJ 

The fixed point scenario suggests a_D + 2aG = 0. A comparison of this con- 
ditions with lattice data is shown in TABLE [21 We observe that the lattice data 
deviate from the theory by about 20%. The Paris and Gent solution [18l [19] 
an = —1, etc = would mean that one needs correction of the vertex renormal- 
ization. The Graz solution [20] ao = — l-2,aG = 0.6 is not compatible with the 
gluon propagator of large lattice [9]. Whether the gluon propagator is affected 
by the two-loop diagram including color antisymmetric ghosts [21] needs to be 
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Table 2. Infrared exponents of MILC and DWF configurations in Landau gauge. 
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Figure 8. The running coupling as{q)/n of MILC/ in Landau gauge The pQCD result of 
Nf = 3 (upper dash-dotted line) and Nf — 2 (lower dashed line) and the extraction of JLab 
are also plotted(blue boxes). 



investigated. I remark that when a large instanton is present, fermions becomes 
massive due to chiral symmetry breaking and they decouple from gluons [43] . 
In this case, the fixed point scenario reflects only an approximate feature of 
IR-QCD, and it could be violated. 

5 Conclusion and discussion 

The color antisymmetric ghost introduces a fluctuation of propagators in the one 
loop level. A comparison of the quenched and the unquenched simulation of SU(3) 
Yang- Mills theory suggests that a quark has the effect of magnifying the square 
norm of the color antisymmetric ghost propagator and reduces its fluctuation. 
The asymptotic behavior of the color diagonal ghost dressing function and 
times the color antisymmetric ghost in the IR is not known. When the crossing 
point of the extrapolation of the color diagonal ghost dressing function and 
times the color antisymmetric ghost is relatively large (~ 60MeV), as in the case 
of DWF, the color diagonal ghost propagator in IR becomes unstable. 

The gluon-ghost-ghost loop contribution with the IR exponent ac = 0.5 
produces the ghost propagator with the dominant component of ac = in 
UV. It is compatible with the recent findings of DSE equation. Historically high 
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Figure 9. The running coupling ai{q)/n of MILC3/ and MILC2/ in Coulomb gauge. The 
pQCD result of Nf = 3 (upper dash-dotted line) and Nf — 2 (lower dashed line) are also 
plotted. 



infrared singularity of the ghost propagator was required to cancel the IR 
singularity of gluon propagator in the tree level [44]. Recently in 3-dimensional 
Landau gauge SU(2) Yang-Mills theory, infrared divergence of the three gluon 
vertex was suggested [38]. It is in conflict with the 4-dimensional Landau gauge 
SU(3) Yang- Mills theory [l6], which asserts that the three gluon vertex is infrared 
vanishing due to instanton effects. If instantons play a role in infrared QCD, 
quenched and unquenched simulation would produce differences, since fermionic 
zero mode which is absent in quenched simulation, is expected to cancel bosonic 
zero mode divergence [47J. 

The IR exponent ac = and = — 1 violates the fixed point scenario 
of the running coupling which requires ao + '^olq = 0. The difference of the 
ghost-gluon coupling of the Coulomb gauge and that of the Landau gauge below 
q ^ 0.6 GeV suggests the vertex renormalization, i.e. the IR suppression of as{q) 
in the Landau gauge may be due to the singularity of the color antisymmetric 
ghost propagator that disturbs the color diagonal part and invalidates the tree 
level approximation. The agreement of the quark- gluon coupling of the Coulomb 
gauge, that of the Landau gauge and the JLab experimental data [U] suggest 
that the different momentum dependence of Landau gauge effective quark-gluon 
coupling and the effective ghost-gluon coupling in the IR is due to our incom- 
plete treatment of the ghost propagator, since I think the universality should be 
preserved in the whole momentum region. 

Since lattice data are taken on the torus and suffers from finite size effects, 
we need to continue careful comparisons of ansatze and results of the DSE and 
the lattice calculation. More extensive study of the Coulomb gauge DWF quark 
propagator |45] of larger lattice is left in the future. 
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